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ABSTRACT 


A  system  that  approximates  the  performance  of  the  optimum  processor 
(optimum  in  the  sense  of  maximizing  array  gain)  is  presented.  Consisting 
of  many  twi. -element  subsystems,  the  system  is  called  a  binary  array  proc¬ 
essor  (BAP).  The  two  main  advantages  of  the  BAP  system  are  that:  (1)  it 
should  be  easier  to  implement  than  the  optimum  processor  because  inver¬ 
sion  of  large  matrices  is  not  required;  and  (2)  It  should  not  suffer  from 
errors  made  in  the  statistical  estimation  of  cross- spectral  densities  of  the 
noise  field  because  the  system  adapts  against  such  errors. These  advantages 
are  obtained  at  the  expense  of  a  few  db  of  array  gain,  but  a  sample  calcu¬ 
lation  illustrates  that  the  array  gain  of  the  BAP  system  can  still  be  signifi¬ 
cantly  greater  than  that  obtained  by  conventional  time-shift-and-sum 
beam  formation. 
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A  SUBOPTIMUM  APPROACH 
TO  ADAPTIVE  ARRAY  PROCESSING 


INTRODUCTION 

The  application  of  optimum  space-time  processing  theory  requires  the 
inversion  of  a  matrix  containing  N2  elements  at  each  frequency,  where  N 
is  the  number  of  hydrophones  in  the  array.  Such  matrix  inversion  introduces 
two  serious  problems:  (1)  The  implementation  of  an  optimum  system  would 
involve  an  impractical  amount  of  hardware;  and  (2)  Since  the  largest  gains 
from  the  optimum  processor  are  obtained  when  the  matrix  approaches  singu¬ 
larity,  any  errors  in  the  statistical  estimates  of  the  matrix  elements  could 
have  a  drastic  effect  on  the  inverse  matrix.  Even  the  inverse  of  an  extremely 
well-conditioned  matrix  could  be  affected  severely  by  estimation  errors  be¬ 
cause  of  the  large  number  of  multiply/add  operations  involved  in  its 
calculation. 


This  report  presents  a  suboptimum  scheme  that  should  not  suffer  from 
measurement  error  because  the  system  adapts  against  estimation  errors  it 
has  made.  Furthermore,  the  largest  matrix  the  system  must  invert  is  a 
2x2  matrix  (such  inversion  is  trivial),  and  hence  the  system  should  be  much 
easier  to  implement  than  the  optimum  processor.  The  only  constraint  on  the 
suboptimum  system  is  that  the  number  of  hydrophones  in  the  array  must  be 
some  power  of  2. 


DISCUSSION 

We  first  assume  that  the  array  has  been  electronically  steered  so  that 
the  desired  signal  is  identical  at  each  input.  In  the  straightforward  approach 
to  optimum  array  processing,  a  linear  filter  is  placed  in  each  channel  (Fig.  1), 
and  then  the  question  is  asked,  "What  must  the  transfer  functions  of  these 
filters  be  in  order  to  maximize  the  array  gain?"  The  required  mathematical 
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Fig.  1  -  Configuration  of  Optimum  Processor 

calculations  then  yield  the  optimum  complex  weights  1  at  each  frequency: 


Z  =  q-  1  a  , 


(1) 


where 

Z  is  an  N  x  1  column  vector  of  the  optimum  weights, 

Q  is  an  N  x  N  matrix  of  the  cross-spectral  densities  between  inputs 
(after  steering),  and 

T 

1 

a  - 

1 

I—  —J 


A  system  that  approximates  the  performance  of  the  optimum  processor  is 
shown  in  Fig.  2.  Consisting  of  many  two-element  subsystems,  the  system  is 
called  a  binary  array  processor  (BAP).  The  BAP  is  constrained  in  that  the 
number  of  hydrophones  must  be  some  power  of  2  ,  but  for  many  applications 
this  constraint  is  not  a  serious  one. 


*D.  J.  Edelblute,  J.  M.  Fisk,  and  G.  L.  Kinnison,  "Criteria  for  Opiimum-Signal-Detection 
Theory  for  Arrays,”  Journal  of  the  Acoustical  Society  of  America,  vol.  41,  no.  1,  January  1967, 
pp.  199-205. 


2 


The  BAP  system  operates  in  the  following  manner:  First  the  subset  of 
filters  {Zj;  f  ,  where  i  =  1,  2,  •  •  •  ,  2M  and  M  =  log2  N,  is  determined  by 
treating  each  two-element  subsystem  separately,  according  to  Eq.  (1).  After 
these  filters  have  been  determined  and  their  outputs  have  reached  a  steady- 
state  condition,  the  subset  {  Z2i },  i  =  1,  2,  •  •  •  ,  2(M‘1),  is  determined  in  the 
same  fashion.  The  procedure  is  repeated  for  the  subset  |z}i  |,  i  =  1,  2,  •  •  •  , 
2<M*2\  and  so  on  until  all  the  filters  have  been  determined.  An  important  fea¬ 
ture  of  the  BAP  system  is  now  evident.  Since  the  filters  |Z2i  }  are  designed 
after  the  filters  j  Zi;| ,  they  adapt  against  errors  in  the  latter  set  caused  by 
incorrect  statistical  estimates  of  cross-spectral  densities.  Similarly,  the  fil¬ 
ters  |  Z  3i}  adapt  against  errors  made  in  determining  the  filters  {  Z2j  | ,  and 
so  on. 


Another  important  feature  of  the  BAP  system  is  that  it  never  requires 
inversion  of  a  large  matrix;  many  2x2  matrices  must  be  inverted,  but  the 
operations  involved  are  trivial.  Hence  the  BAP  system  should  be  orders  of 
magnitude  simpler  to  implement  than  the  optimum  system. 


ANALYSIS  OF  THE  BAP  SYSTEM 


Although  it  is  possible  to  analyze  the  BAP  system  for  the  general  case,  we 
car.  simplify  the  analysis  greatly  by  assuming  that  the  input-noise  cross- 
spectral  density  matrix,  in  addition  to  being  Hermitian,  is  a  Toeplitz  matrix, 
that  is,  it  has  equal  elements  along  any  diagonal.  Only  this  special  case,  for 
which  the  transfer  functions  of  the  filters  in  each  binary  subsystem  are  complex 
conjugates,  is  examined  in  this  report.  Many  filters  have  identical  transfer 
functions,  as  given  by  the  following  relation: 

Zki-Zk(U2>-  k  =  1,  2,  •  •  • ,  (M  —  1  ] 

i  =■  1,  2,  •  •  • ,  [2(M-k  ♦!)  _  2]  . 

Figure  3  shows  the  BAP  processor  configuration  for  the  special  case  under 
consideration. 
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Fig.  3  -  Binary  Array  Processor  for  Case  when  Input  Cross-Spectral  Density  Matrix  Is  Toeplitz 


Before  proceeding  further  with  the  analysis,  we  must  present  the  notation 
that  will  be  used  for  all  cross-spectral  density  matrices.  At  the  input,  i.  e. , 
just  before  the  filters  j  Z,j| ,  we  have: 


QA  = 


QAn 


Qa1(2m) 


QA, 


(2M),  •  ’  ’  QA{2M)(2M) 


and  just  before  the  filters  |  Z2i  |  , 
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QB  = 


QB 


1 1 


QB  if  ;M- 


l(jM-l) 


The  matrix  of  spectra  preceding  the  filters  j  Z  3i|  is  QC,  and  so  on. 

Analysis  of  the  BAP  system  proceeds  as  follows.  Given  QA,  the  filters 
jZ(i}  are  determined,  and  then  the  matrix  QB  is  found;  similarly  the  fil¬ 
ters  |Z2i|  and  the  matrix  QC  are  determined.  The  process  continues  until 
the  power  spectrum  of  the  final  output  stage  is  known.  The  array  gain  can  then 
be  c&lculatec  after  this  same  procedure  is  repeated  using  the  same  filters  with 
signal  (instead  of  noise)  cross-spectral  density  matrices. 

When  QA  is  a  Hermitian,  Toeplitz  matrix,  the  filters  {zu|  are  all  deter¬ 
mined  once  ZM  is  known,  and  Z;1  is  found  by  solving  the  following  matrix 
equation: 


- 

- 

-i 

■  * 

7,, 

QAn 

i 

• 

Zn 

1 - 

o 

> 

_  • 

< 

O' 

i 

. . 

Now  since  QB  is  also  a  Hermitian,  Toeplitz  matrix,  it  is  necessary  only 
to  find  QB QB12,  QB1},  •  •  •  ,  QB,(2m-i>  in  order  to  determine  QB.  The 
element  QBn  is  given  b\  the  following  expression  from  linear  system  theory *: 


QB , ,  =  QA , ,  |  Z 


1 1 


2  +  QA,2Z*i  +  QA21zj,  .QA22  | Z 


I  1 


Jj.  Bcndat  and  A.  Piersol,  Measurement  and  Analyst!  of  Random  Data,  John  Wiley  and  Sons, 
N.  Y.,  1966,  Eq.  9.166,  p.  108. 
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With  the  use  of  the  relations 


QA,,  -  QA2, 
and 

0A2,-0A12 


we  reduce  this  equation  to 

.  |  ■-  ’tM ,  |  2  ♦  2Rc  jgAu/’j 


(2) 


The  element  QBJ2  is  derived  here  because  no  expression  is  readily  avail¬ 
able  in  the  literature.  In  Fig.  4  the  filters  of  the  first  four  channels  are 


NA,  (f ) 

na2(») 
na3  (t) 
na4  ( r) 


Fig.  4  -  Derivation  of  QB  12 


represented  by  their  impulse  responses  instead  of  their  transfer  functions. 
The  functions  NA,(t),  NA2(t),  NA  ,  (t),  NA  4(t),  N  B ,  (t) ,  and  NB2<t)are 
noise  timt  functions.  The  expressions  for  NB,(t)  and  NB2(t)  are 


and 


NB  ,(t)  = 


NB  2(t)  = 


(z  J  (  (at  NA  j(t  -  o)  ♦  z  ( ,(<?) N  A  2(t  -a)]  do 


(z  j  ,  (p)  NA  ,(i  -  p)  ♦  z  l:(p)NA4(l  -  p)l  dp  , 
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and  the  cross-correlation  function  between  these  two  functions  is 


RB[2(t)  =  NB,(t)NB2(t  ♦  t)  , 


where  the  bar  denotes  statistical  averaging.  Substitution  of  the  expressions  for 
NB^t)  and  NBJt)  into  this  equation  gives 

<v> 

RB(2(7)  (zjjialNAjit-fj)  +  Zp(a)NA,(t -a)llzn(p)NA3(t -p  +  r)  +  z^iplNA^it  -  t)1  dadp 


-no  -  Qe 


RB 


ne  OC 

,( T )  =  I  I  [zjjialzjjiplRApiT  +  o  -  p)  *  z12(o)zn(p)RA  2}(~  +  a-p) 

•ou  -Ou 

+  z  ||  1c)  z  jj(p)  R  (4  (T  +  o  -  p)  +  zn(o)zl2(p)  RA  24  (T  +  a  -  p)l  da  dp  , 


where 


RA  m(t+  ct-  p)  =  NA,(t  -  a)NA,(t  -  p  +  r)  , 


RA^fTt  a  -  p)  =  NA2(t-CT)NA}(t-p  +  t)  , 


RA  m(t  ♦  a  -  p)  =  NA  j(t  -  ct)NA4(i  -  p  +  t)  ,  and 


RA24(t  +  ct  -  p)  =  NA2(t  -  ct)  N  A  4(t  -  p  +  T)  . 


The  Fourier  transform  of  RBl2(r)  is  the  quantity  we  desire: 

QB  ,  2  =  (  RB12(T)e-,JnfT  d  T 


Substituting  the  previously  given  expression  for  RB12(t)  into  the  above  integral 
and  multiplying  and  dividing  the  integrand  by  e ->jT,f(CT-p)  yields 


no  no  no 


QB 


1  2 


j"  j  |  ^uMcnv(azJ](p)e  '2n(pRAli(T  +  o-p)e-'2n(lr*<T-f,)  dodpdi 


-no  -no  -no 


z  1 2<a)  c 


i  2”  Icr 


1 1 


(p)  c"1  2tt^  RA 


2} 


(r  ♦  a  -  p) 


•j  2nf(rt  a  -  p) 


dadpdr 


no  no  no* 

r  i 

-no  -no  -no 


z..(o)''2vic’z.Ap)'-'2n,cRA.\7  +  a-p) 


-\l-nRr*o-p) 


dadpdr 


no  no  no 


iff  z12(a)ciJT’fCTz12(p)t'2”f',RA24(r  >  a  -  p)  e’i  mHr*a-p)  d<7(JpdT  , 


-no  -no  -a 


which  reduces  to 


QB,2  =  7.;i(f)Zn(f)QA1,  4  Z*2(f)7.n(f)QA2, 

4  /.;i(f)ZI2(f)QAM  4  Z*2(f)Z12(f)OA24  ■ 


Now  since 


Z,2(f) 

z,i<f), 

<M;, 

Qa,2, 

and 

<^24 

yA13, 

the  element  QB[2  becomes 

QB,2  2<Mn 

1  ^  1  1  1  2 

♦  QAjjZ] [  4  QA)4Zi i  . 

The  expressions  for  QBU  , 

QBm 

,  etc.  are  obtained  from  the  following 

general  expression,  which  was  derived  by  changing  subscripts  in  the  preceding 
derivation: 
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kM . 


(3) 


-  *2 

QH,K  2QA1(2K-1)  Z  I  j  *  QA,(2k_2)/.i  ,  +  QA  1(2K)^1  !  • 

Thus  all  the  elements  of  the  matrix  QB  can  be  determined.  In  fact,  we  can 
now  determine  the  entire  BAP  system  simply  by  changing  subscripts  and  varia¬ 
ble  names.  For  example,  the  filters  jz2i|  are  found  in  the  same  manner  as 
the  filters  |Zi;| ,  the  only  difference  being  that  elements  of  QB  are  used  in¬ 
stead  of  elements  of  QA.  Similarly,  the  elements  of  QC  are  obtained  from 
the  same  expressions  used  for  QB  simply  by  changing  all  the  A's  to  B's  and 
all  the  B's  to  C's.  These  procedures  are  repeated  until  all  the  filters  are 
determined  and  the  value  of  the  output-noise  power  spectral  density  is  found. 
Then  the  same  filters  and  equations  with  noise  cross-spectra  replaced  by 
signal  cross-spectra  are  used  to  find  the  value  of  the  output-signal  power  spec¬ 
tral  density.  Array  gain  is  then  computed  by  dividing  the  ratio  of  signal  power 
spectral  density  to  noise  power  spectral  density  at  the  output  by  the  same  ratio 
at  an  input  reference  hydrophone. 


SAMPLE  CALCULATION 

Consider  a  linear  array  oriented  vertically  in  a  surface-generated  noise 
field  and  having  four  elements.  The  following  expression,3  which  has  been 
normalized  with  respect  to  a  selected  input  reference  quantity  QAU  .  defines 
the  elements  of  the  matrix  QA  : 

sin  R  cos  fl  „ 

r  mo  f  mn 

In  this  expression. 


d_„  is  the  distance  between  the  mth  and  nth  receivers  in  feet,  and 
mn 

A  is  the  wavelength  in  feet  . 


’B.  F.  Cron  and  C.  H.  Sherman,  "Addendum:  Spatial-Correlation  Functions  for  Various 
Noise  Models,"  Journal  of  the  Acoustical  Society  of  America,  vol.  38,  no.  5,  November  1965, 
p.  885. 
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Since  we  want  QA  to  be  Toeplitz,  let  d^/A  be  constant  for  all  (m,n).  For 
dmn/A  =  0.3,  the  elements  that  specif\  the  matrix  QA  are: 


The  matrix  QB  is  determined  once  QBn  and  QB,2  are  known,  and  they  are 
found  from  Eqs.  (2)  and  (3); 


<?nn 

2QA,, 

2Rc  [QA12  /.*2  1  -  0.20 

and 

QB,2 

2QAI3 

Z 1 1  2< 

QA  i  j  /. ,  f  +  j  4  Z  n 

QB,2 

0.025  - 

j(0. 155) 

0.155  /-80. 85°  . 

The  filters  jZ2i|  are  all  found  from  Z2J,  which  is  given  by 


Z2i  =  QBn  -  QB i j  0.235  ♦  j(0. 155)  =  0.282  /33.4°  , 

and  the  noise-output  power  spectrum  is  found  from  Eq.  (2)  to  be 

QC  2QBn  jz21i2  ♦  2Rc  [QB^Zjjj  =  0.203  . 

The  above  procedure  is  now  repeated  with  the  same  filters  but  with  signal 
cross-spectra  in  place  of  the  noise  cross-spectra.  We  have  assumed  that  the 
signal  is  identical  at  each  input;  therefore,  the  input-signal  cross-spectral 
density  matrix,  normalized  with  respect  to  SAU  ,  is 

1111 
1111 
1111 
1111 
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The  elements  SBn  and  SB12  are  found  from  Eqs.  (2)  and  (3)  with  SAn, 

SA12,  SA,j,  and  SAU  replacing  QAUl  QA12,  QA15,  and  QAm,  respec¬ 
tively.  They  are 

SB ,  j  2SAn  Zn|2+  2Re  |SA|2Z*  J  2. 14A 
and 

SB12  2SA,,  Z,,  *\  SA12/.-;,  .  SA,/.;;  2.143. 

Thus  the  matrix  SB  is  determined.  Equation  (2)  with  the  appropriate  substitu¬ 
tions  then  yields  the  signal-output  power  spectrum  as: 

SCn  2SB,,  1Z21  l2  ♦  2  Re  (SB,//*  J  4.74. 

The  array  gain  of  the  system  is,  therefore, 

sc 

G  =  1 — LL  23.1  =  13.6  db  . 

O'  n 

The  above  sample  calculation  has  showed  an  array  gain  of  13.6  db  for  a 
four-element  BAP  system  operating  with  a  vertical  array  in  a  surface-generated 
noise  field.  Optimum  processing  under  the  same  conditions  would  give  an  array 
gain  of  19  db,  whereas  conventional  time-shift-and-sum  beam  formation  would 
give  less  then  10  log10(4)  =  6  db.  Thus,  although  the  BAP  system  may  not  re¬ 
sult  in  as  large  an  array  gain  as  the  optimum  processor,  it  still  provides  sig¬ 
nificantly  higher  gains  than  are  obtainable  by  conventional  methods. 


CONCLUSIONS 

A  suboptimum  approach  to  array  processing  has  been  described.  The  sys¬ 
tem,  called  a  "binary  array  processor"  (BAP),  consists  of  a  collection  of  many 
two-element  subsystems  The  relative  ease  of  BAP  implementation  (resulting 
from  the  fact  that  inversion  of  large  matrices  is  not  required)  should  compen¬ 
sate  for  the  sacrifice  of  array  gain.  In  addition,  the  BAP  system  should  be  less 
sensitive  than  the  optimum  processor  to  estimation  errors  in  the  noise  cross- 
spectral  densities  because  it  adapts  against  such  errors. 
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essor  (BAP).  The  two  main  advantages  of  the  BAP  system  are  that:  (1)  It 
should  be  easier  to  implement  than  the  optimum  processor  because  inversion 
of  large  matrices  is  not  required;  and  (2)  It  should  not  suffer  from  errors 
made  in  the  statistical  estimation  of  cross -spectral  densities  of  the  noise 
field  because  the  system  adapts  against  such  errors.  These  advantages  are 
obtained  at  the  expense  of  a  few  db  of  array  gain,  but  a  sample  calculation 
illustrates  that  the  array  gain  of  the  BAP  system  can  still  be  significantly 
greater  than  that  obtained  by  conventional  time-shift-and-sum  beam  formation. 
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